Edwards entropy and compactivity in a model of granular matter. 
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Formulating a statistical mechanics for granular matter remains a significant challenge, in part, 
due to the difficulty associated with a complete characterization of the systems under study. We 
present a fully characterized model of a granular material consisting of N two-dimensional, fric- 
tionless, hard discs, confined between hard walls, including a complete enumeration of all possible 
jammed structures. We show the properties of the jammed packings are independent of the dis- 
tribution of defects within the system and that all the packings are isostatic. This suggests the 
assumption of equal probability for states of equal volume, which provides one possible way of 
constructing the equivalent of a microcanonical ensemble, is likely to be vaild for our model. An 
application of the second law of thermodynamics involving two subsystems in contact shows that 
the expected spontaneous equilibration of defects between the two is accompanied by an increase in 
entropy and that the equilibirum, obtained by entropy maximization, is characterized by the equal- 
ity of compactivities. Finally, we explore the properties of the equivalent to the canonical ensemble 
for this system. 

PACS numbers: 81.05.Rm, 83.80.Fg, 05.20.Jj 
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I. INTRODUCTION 

Granular materials such as sand are too massive to be 
influenced by the thermal fluctuations that affect the mo- 
tion of particles on the atomic scale. If left undisturbed, 
a sand pile remains in a rigid, solid-like structure, but 
when poured, sand flows just like a fluid. Foams, made 
up of a collection of macroscopic bubbles, can support a 
yield stress like a solid but they flow like a fluid when 
sheared above a certain threshold rate. In both cases, 
the particle-particle structure and particle dynamics are 
similar to that observed in atomic and molecular liquids 
and glasses, suggesting that many athermal systems may 
exhibit a wide variety of physical phenomena that were 
originally thought to occur only in thermal systems [l[. 
While there have been a number of efforts to capture the 
connection between the different systems, including the 
formulation of a phase diagram for jammed and glassy 
materials 0, we do not have an underlying and gen- 
eral theory to describe granular matter. The machinery 
of classical thermodynamics and thermal statistical me- 
chanics, which provides us with the tools for the study of 
complex behavior in molecular systems, cannot be used 
to explain the behavior of a granular system because 
the energy required to move a macroscopic particle far 
exceeds that available from thermal excitation. Conse- 
quently, there is an extensive effort to develop a form of 
statistical mechanics that can account for the properties 
of athermal systems . 

An appropriate starting point for a statistical mechan- 
ics of granular systems is to develop a probability distri- 
bution function that is equivalent to the microcanonical 
ensemble where all the states of the same energy are as- 
sumed to have an equal probability of being sampled. 
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One possibility, suggested by Edwards 0, is to as- 
sume that all jammed states of equal volume have equal 
probabilities so that the probability of finding a given 
configuration i is, 



P t = e~ s ' x 8{V - Wi)Q, 



(1) 



where V is the volume, Wi is the Hamiltonian-like vol- 
ume function for the configuration, S is the entropy and 
A is the analogue of the Boltzmann constant. The nor- 
malization in Eq. [T] is 



s/x 



S(V - Wi)Qd(&U degrees of freedom), (2) 



and 
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1 if collectively jammed 
otherwise 



(3) 



Evaluating the integral in Eq. [2] and taking the natural 
log of both sides defines the entropy as, 



s/x = inn{v), 



(4) 



where ^(VQ is the number of collectively jammed configu- 
rations [10l4l2| with a given V. il(V) grows exponentially 
with the number of particles, N, so that entropy is ex- 
tensive and the compactivity, which is equivalent to the 
temperature in a thermal system, can be defined as 



X 



dV 

as' 



(5) 



In general, fi(V) is not known. One difficulty is 
that little is known about how particles share volume in 
jammed materials and while a number of different tessel- 
lation methods have been suggested [&-:8j, the construc- 
tion of W remains a significant challenge. Furthermore, 
the global nature of the collective jamming condition 
Il2j does not allow us to determine which configurations 
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are jammed on the basis of the local geometric proper- 
ties of the individual particles. The effects of friction and 
stress in jammed configurations add further complica- 
tions, leading to the proposal of alternative distribution 
functions @. 

A numerical simulation of a slowly, sheared granular 
material [l3T ] suggests a thermodynamic temperature, re- 
lated to compactivity, can be extracted from particle mo- 
bility measurements, but simulations of frictionless hard 
discs 14j appear to cast doubt on the assumption that 
packings o f eq ual volume are equi-probable. Recent ex- 
periments [l5| on externally-agitated granular systems 
found that while compactivity did not equilibrate be- 
tween two subsystems, both subsystems shared a num- 
ber of reproducible properties, strongly supporting the 
notion that a thermodynamic approach to granular ma- 
terials is a real possibility, even though we do not have a 
complete picture of how this should be developed. 

However, for all the systems studied to-date, the dis- 
tribution function and a detailed knowledge of the par- 
ticle packings is lacking. The goal of the present work 
is to examine the statistical mechanics of a simple model 
granular system that has the potential to be examined ex- 
perimentally and for which the full density states can be 
obtained exactly. We analyze the distribution of jammed 
packings for a system of highly confined hard discs in Sec- 
tion II and show that the assumption of equal probabil- 
ity for states of equal volume is likely to be valid for this 
model. In Section III, we examine the predictions of the 
second law of thermodynamics by considering the equi- 
librium state obtained by bringing together two isolated 
systems into the granular equivalent of thermal contact. 
The canonical ensemble is studied in Section IV while 
Section V contains our discussion and conclusions. 



particles are only able to contact their nearest neighbor 
on each side and the wall so there are only two local pack- 
ing environments. In addition, the confinement prevents 
collective motions that would allow the locally jammed 
structures to collapse and the complete distribution of 
collectively jammed states can be obtained from local 
geometric considerations alone. 

Fig. QJ, shows a jammed configuration containing the 
two possible local disc arrangements. The most dense 
state involves two particles contacting across the chan- 
nel while the defect state consists of two neighboring 
particles contacting along the channel, which results in 
a less dense arrangement. Both the most dense local 
arrangement and the defect have well defined volumes, 
vq = H^J (2cr — H)H and v\ = Ha respectively, that are 
additive so the volume function, W, can be written as, 



V = W = (N - M)v + Mvx 
= NH \y/(2 ~ H)H(1 



(6) 



where M is the number of defects in the system, 9 = 
AI/N is the fraction of defects and a is the unit of 
length. The occupied volume fraction for a packing is 
4>j = ir/(4H[y/(2-H)H(l - 9) + 9}). _ 

To count the number of packings with a given V, we 
develop a lattice gas, or Ising model-like, description of 
each configuration of discs by drawing a bond between 
the centers of neighboring discs and labelling the defect 
bonds as "1" and the most dense bonds as "0" (see Fig[l}. 
This means the total number of bonds equals N. Follow- 
ing Hill if M is the number of defect bonds, we can 
divide a configuration into blocks consisting of consecu- 
tive ones in the chain and blocks consisting of consecutive 
zeros in the chain. A block of ones is necessarily sepa- 



II. PACKINGS OF HIGHLY CONFINED HARD 
DISCS 

Our model granular material consists of N two- 
dimensional (2D) frictionless, hard discs, with diameter 
a, trapped between two hard walls separated by a dis- 
tance H/a < 1 + \/3/4 and we will ignore the effects of 
gravity. The density of states for this system was previ- 
ously described in Ref. [l6[ and is included here for the 
sake of completeness and to highlight the fact that the 
distribution of packings for this system can be obtained 
using simple combinatorial arguments. In addition, we 
also show here that all the packings of the system are 
isostatic and that packings of the same volume are all 
structurally equivalent. 

In 2D, a particle is locally jammed if it has at least 
three rigid contacts that are not all in the same semi- 
circle. However, while local jamming of all particles is a 
necessary requirement for a collectively jammed state, it 
is not sufficient because the concerted motion of a group 
of particles can cause a packing to collapse. By confining 
the discs between walls separated by H/a < 1 + \/3/4, 
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FIG. 1. (a) A jammed configuration. The dashed lines join 
the centers of contacting discs and are labelled 0, for the most 
dense, or 1 for the least dense (defect) arrangement. The 
shaded regions represent the volumes associated with each, 
(b) An arrangement with two neighboring defects (-1-1-) is 
not jammed because the shaded disc does not fulfill the local 
jamming condition and is able to move vertically. 
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rated from a block of zero by a boundary consisting of a 
1-0 or 0-1 bond. The total number of configurations with 
M defects and M i boundaries is obtained by considering 
the number of different ways of arranging M ones among 
the (Moi + l)/2 possible block of ones, such that there 
is at least one "1" in a block, and (N — M) zeros among 
the (Mqi + l)/2 possible block of zeros. This gives, 



fij(M,M i) 



M\(N — M)\ 



[M- ^fi ]![iV- M - ^p-]![^f!] 2; 



(7) 



in the limit of large numbers. The volume of a packing of 
N discs only depends on = M/N and is independent of 
Moi. If the defects could be distributed totally randomly 
throughout the system, then Eq. [7] would reduce to the 
usual Ising model expression, flj(V) = N\/M\(N — M)\. 
However, not all defect arrangements result in truly 
jammed states. When two defects appear next to each 
other in a "1-1" arrangement, the central particle does 
not satisfy the local jamming condition (Fig. \\)p). This 
arrangement represents a saddle point in configuration 
space that will collapse if the central particle is per- 
turbed in a direction normal to the wall. These con- 
figurations can be eliminated from the distribution by 
setting Moi = 2M, which ensures that every defect is 
isolated from the other defects and gives, 



(N-M)\ 
M\(N-2M)V 



(8) 



The distribution of packings is binomial with a single 
most dense structure, containing no defects, a single least 
dense structure with = 0.5 and a maximum number of 
packings when 9 — 1/2 — a/5/10. 

Mechanical equilibrium is an important property of 
granular packings that can be understood in terms of 
a balance between the total degrees of freedom in a sys- 
tem and the number of force equations that constrain or 
counteract them [Tj|. In our system, each disc has two 
translational degrees of freedom, resulting in a total of 
2N degrees of freedom. Each disc in a jammed packing 
contacts two other discs and the wall, which results in a 
total of N disc - disc force equations, noting each con- 
tact involves two discs, and N disc - wall force equations, 
giving a total of 2N equations. The number of degrees 
of freedom is exactly balanced by the number of force 
equations which shows that all the jammed packings for 
this model are isostatic, independent of the number of 
defects. Furthermore, since each particle only contacts 
one of its neighbors on either side, it does not feel any 
interaction from particles beyond its nearest neighbours 
so there is no interaction between defect states. As a re- 
sult, the properties of the collectively jammed packings, 
with the same volume, should be independent of how the 
defects are distributed i.e. there is no distinction between 
a packing with an ordered arrangement of defects com- 
pared to one with a more random arrangement. These 
properties strongly suggest that the assumption of equal 
probability for states of equal volume is valid for the cur- 
rent model. 



III. THE SECOND LAW OF 
THERMODYNAMICS 

The second law, in classical thermodynamics, postu- 
lates the existence of entropy as a state function in order 
to understand the driving force behind spontaneous pro- 
cesses. Furthermore, it states that dS > for any spon- 
taneous process in an isolated system so that the entropy 
is maximized when the system is in equilibrium. The 
most widely used example of the 2 nd law involves bring- 
ing two metal bars at different initial temperatures, say 
T\ <T 2: into thermal contact. If the composite system is 
isolated, then experience (experiment) tells us that heat 
will flow from the hot bar into the cold bar until it comes 
to equilibrium and the temperature in each bar is the 
same. 

We can construct an analogous experiment with our 
granular system as follows: First we note, on the basis of 
Eqs.0]and[51 that the entropy of the system with a given 
volume is, 



iln(l - 9) - (1 - 20)ln(l - 29) 



-01n0, 



and the compactivity is, 



(9) 



X 



dV 
~dS 



dVd9_ 
~d9dS 



HU{2-H)H-1 



ln(l 



ln(0) - 21n(l - 20)' 



(10) 



where we have made use of the derivatives of Eqs. [6] and 
|9]with respect to 0. In the lower limit, X — > as — > 0, 
i.e. when the system is in its most dense state. In the 
upper limit, X — > oo as 1/2 — a/5/10, which occurs 
at the maximum in the distribution of jammed packings. 
Packings with higher concentrations of defects are not 
sampled in equilibrium because these would give rise to 
negative compactivities. 

To prepare two subsystems of N% and N 2 discs, at com- 
pactivities X\ and X 2 respectively, we must in principle, 
place each subsystem in contact with a "thermal" reser- 
voir with which it can exchange volume until equilibrium 
is reached. However, since both the volume and com- 
pactivity are functions of the fraction of defects, it is 
sufficient to select starting configurations with the ini- 
tial defect concentrations 9[ and 6' 2 , which correspond 
to some initial conditions with X\ ^ X 2 . We can now 
bring the two subsystems into contact so that they can 
exchange volume. This can be achieved by placing the 
last particle of system one and the first particle of system 
two in contact so that the configuration remains jammed. 
In the context of the present experiment, an isolated sys- 
tem is one that cannot exchange volume with an external 
reservoir so the composite system has a fixed total vol- 
ume Vt = Vi + V2, where we use large enough system 
sizes that the small volumes associated with the end ef- 
fects can be ignored. 
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The challenge of dealing with a realistic set of dynam- 
ics, such as shaking, will be address in more detail in 
our discussion, but in the meantime we need to use an 
idealized dynamics that will allow the system to jump 
from jammed state to jammed state under the condi- 
tions of fixed volume. As the system is perturbed, we 
expect it to move between jammed states by a series of 
random local particle rearrangements which will require 
the movement of defects. If all the particles are made 
from the same material, our physical intuition tells us 
that the defects will spontaneously move from a region 
of high concentration to a region of low concentration, 
until eventually the system comes to equilibrium and the 
defects are, on average, equally distributed throughout. 
Since both subsystems are made from the same material 
and have the same size, a uniform distribution of defects 
also corresponds to a uniform compactivity. 

Application of the second law of thermodynamics im- 
plies that the equilibrium distribution of defects can be 
obtained by maximizing the total entropy of the compos- 
ite system, St = Si(#i) + 82(02)- The conservation of 
volume gives 



JVi(01-0i 



N 2 8' 2 



No 



(11) 



which can be used in Eq. [9] to yield an expression for 
St solely in terms of 8\ and the initial defect concentra- 
tions consistent with the starting compactivities for each 
subsystem. The equilibrium value of 0\ in the composite 
system, obtained from oISt/oIOi = 0, is then 



n = 



N 2 e' 



Nx +N 2 



(12) 



This is just the average number of defects, or uniform dis- 
tribution, which is consistent with our physical expecta- 
tion. Furthermore, the compactivity in both subsystems 
are the same at equilibrium. 

The present model can also be used to study the ef- 
fects of bringing different types of granular materials 
into contact. The simplest example to study is the 
case where a subsystem of Ni discs with diameter o~\ 
is brought into contact with a second subsystem of N2 
discs with diameter 02, under the conditions ai < 02 
and H/ai < 1 + \/3/4, which ensures the interface be- 
tween the two subsystems can still be made through a 
single nearest neighbor contact. V\ is given by Eq. [6] 
with N = N\ and the volume of subsection 2 is, 



V 2 = N 2 H 



v/(2a 2 - H)H(l 



<J20 



At a fixed Vr, we now have, 



Ni (0i - 9[) y/{2 - H)H - 1 



N 



(V (2<r 2 - H)H - 0-2 



(14) 



which reduces to Eq. [TT] as er 2 °~i = 1 an d can be 
used in Eq. [5] to provide an expression for 52 in terms 
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FIG. 2. (top) The volume per particle for each subsystem, 
V1/N1 (dash-dot line), V2/N2 (dashed line) and the composite 
system, Vt/N (dash-dot line), (bottom) The entropy per 
particle, S1/N1 (dash-dot line) S2/N2 (dashed line) and the 
composite system, St/N (dash-dot line). The vertical dotted 
line locates the maximum in St, Ni = N2 — 100, o\ = 1, 
a 2 = l.lo-i, H = 1.8<7i, 9[ = and 6' 2 = 0.27 



of 0\ . Fig. [5] shows the volume and the entropy per par- 
ticle for the two subsystems and the composite system 
with the initial conditions chosen so that 0[ = and 
0' 2 = 0.27 which corresponds to X\ — and X 2 = 00, 
respectively. At equilibrium, where St is at a maximum, 
X\ = X2 = 0.318, but neither the volumes or entropies 
of the two subsystem are equal because the change in 
volume associated a defect is different in each subsystem 
ms a result of the different particle sizes. 

To show that the compactivities in the two subsystems 
are equal when St is maximized at fixed Vr, we begin 
by assuming X\ ^ X2 and show that this leads a con- 
tradiction to the conservation of the total volume. From 
Eq. [TU] we have, 

_ fdvA fdo 1 \ (dv 2 \ (do 2 \ fd0A 

(15) 

where we have used the fact that 02 is a function of 9\ at 
fixed Vt- Maximizing the total entropy yields 



dSx 
d0i 



dS 2 
80i 



(16) 



(13) 

which, in combination with Eq. [T5l gives. 



dVA _ fdV2_\ (802 
d8 1 r \d02 \d0i 



90i 



^ ■ (17) 



This completes the proof as this last inequality suggests 
that the volume lost from one system is not equal to 
the volume gained by the other as the number of defects 
in system 1 is varied which violates the conservation of 
volume. Therefore X\ must equal X2 at equilibrium. 
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FIG. 3. (a) The entropy, S, and (b) the volume relative to 
the volume of the most dense state, V/V m in, as a function 
X for a system with particles H/a — 1.86 (solid line) and 
H/a = 1.50 (dashed line). 



IV. CANONICAL EQUILIBRIUM 

The equivalent to the canonical ensemble in a thermal 
system is obtained by integrating over all possible degrees 
of freedom and providing a Boltzmann type weight to 
each jammed state on the basis of its volume, which for 
the present system gives, 

= J <de- w/xx d(&\\ degrees of freedom) 

M=N/2 

= tt(M,N;V)e- w/xx 



M=0 

e *x 2-F1 



1 N N AT A 

— — — ■ — A/ • — dp xx xx 

2 2' 2' ' 



a») 



where Y is the analog to the free energy, 2-F1 is a hyper- 
genometric function and the integral has been reduced to 
a sum over the discrete set of jammed states. The free 
energy can now be used in conjunction with the equiva- 
lent thermodynamic machinery to obtain other thermo- 
dynamic variables. For example, in analogy with the 
Helmholtz free energy, given 

Y = V - XS, (19) 

the entropy of the system can obtained from 



Fig. [3] shows the entropy and volume relative to the 
most dense packing for a granular system of discs con- 
fined within channels of diameter H/a = 1.86 and H/a = 
1.50 as a function of the compactivity. The underly- 
ing distributions of jammed states is the same for both 
systems and the difference in behavior between the two 
arises from the larger change in volume associated with 
the gain or loss of a defect state in the wider channel. 
As a result, the relative change in volume is greater, and 
occurs gradually as a function of X, for the wider chan- 
nel. The properties of the narrow channel are relatively 
insensitive to changes in X at the higher compactivities, 
but we see a rapid decrease in entropy and volume as X 
approaches zero. 



V. DISCUSSION 

The goal of the present paper was to explore the statis- 
tical mechanics of a simple model granular system. We 
are able to carry out a complete analysis of the distribu- 
tion of jammed states for a system of highly confined hard 
discs and find that all the packings of a fixed volume are 
structurally the same. This suggests that the Edwards' 
assumption of equi-probability for jammed states of equal 
volume is well justified in the present model, even if it is 
not valid for more complex granular systems. The model 
also represents a genuine granular system that can be re- 
alized experimentally and it is rare to find cases where 
the exact results for real experiments are available. 

However, there arc aspects of granular systems not cov- 
ered in this model. First, we point out that our analysis 
required an idealistic dynamics that allowed the system 
to move directly between jammed states at fixed volume. 
Shaken granular materials need to generate a momentary 
increase in free volume to allow the particles to unjam, 
move and then return to a jammed state. The role of the 
energy injection and how the free volume affects the way 
the set of jammed states is sample is a critical feature of 
granular systems that has not been included here. The 
presence of friction between particles has also been ig- 
nored and we would expect this to alter the number and 
type of packings included in the distribution, as well as 
possibly influencing the assumption of equi-probable vol- 
ume states. For example, it is immediately obvious that 
the neigboring defect state (Fig.[T|3) would become stable 
in the presence of friction, but other "bridge-type" parti- 
cle contacts may also need to be included. Nevertheless, 
the simplicity of the present system and its accessibilty 
to experiment, suggest that it may be an invaluable tool 
in understanding some of the elementary aspects in the 
statistical mechanics of granular systems. 
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